KA-TP-1-1996 
January 1996 



The Measure of the Orthogonal Polynomials Related to 
Fibonacci Chains: The Periodic Case * 



Wolfdieter Lang ** 

Institut fiir Theoretische Physik 
Universitdt Karlsruhe 
Kaiserstrasse 12, D-76128 Karlsruhe 
Germany 



ABSTRACT 

The spectral measure for the two famihes of orthogonal polynomial systems related to 
periodic chains with N-particle elementary unit and nearest neighbour harmonic interaction 
is computed using two different methods. The interest is in the orthogonal polynomials 
related to Fibonacci chains in the periodic approximation. The relation of the measure to 
appropriately defined Green's functions is established. 

PACS numbers: 02.10Nj, 33.10Gx, 61.44.+p 



* Herrn Dr. Harald Fried in Dankbarkeit gewidmet 

** E-mail address: wolfdieter.lang@phys.uni-karlsruhe.de 



2 



1. Introduction 

Two associated families of orthogonal polynomial systems govern longitudinal time- 
stationary vibrations of linear chains with nearest- neighbour harmonic interaction (springs 
Kn, masses rrin)- 

These polynomials are defined by three-term recurrence relations appropriate for or- 
thogonal polynomials [1]. In the case of a mono-atomic chain (mass mo) with uniform 
coupling (strenght k) they are deformations of Chebyshev's S'^(2(l — x)) = Un{l — x) 
polynomials of the second kind, x = u!'^/2u>q, with ujq = K/mo, is a normalized frequency- 
squared. 

Our interest is in Fibonacci-chains which have uniform coupling («;„ = k) and two 
masses (mass- ratio r = mi /mo) distributed at site number n = 1,2,... in accordance 
with the binary sequence 1,0, 1,1,0, 1,0,... This quasi-periodic sequence is generated by 
the Fibonacci substitution rule 1 1,0 and 0^1. Such chains have been considered 
as simple models for special binary alloys [2]. The same structure is encountered in the 
problem of the phonon spectrum of a one-dimensional Fibonacci quasi crystal [3]. In this 

case the associated orthogonal polynomial systems are denoted by {Sn\x)} and {Sn\x)} 
[4] . Up to now their spectral measures (or moment functionals) have not been determined. 
From rigorous results on the Fibonacci Hamiltonian in the context of the one- dimensional 
Schrodinger equation one expects that this measure is of the singular continuous type 
supported by a Cantor set of zero Lebesgue measure [5]. 

As an approximation to the quasiperiodic problem we identify in this work the spectral 
measure for periodic Fibonacci chains with an elementary unit consisting of N masses 
following the pattern of the first N entries of the above given binary sequence. The 
measures for general A^-periodic orthogonal polynomial systems {Sn{x)} and {Sn{x)}, 
defined by three-term recursion relations with A— periodic coefficients, can be found using 
two different methods. 

The first method employs the Bloch-Floquet solutions for these periodic chains and is 
based on the evaluation of a (judiciously chosen) complex contour integral. This is a special 
application of a general method available for orthogonal polynomials with asymptotically 
periodic coefficients in the recursion relation [6] . For our purpose only the strictly periodic 
case is of interest. With this restriction a similar procedure is discussed in ref.[7]. The 
second method uses the fact that the continued fraction associated with orthogonal poly- 
nomials is the Stieltjes transform [8] of the spectral measure [lb, 9]. For periodic systems 
this continued fraction can be determined as fixed-point of a certain Mobius transfor- 
mation, and the measure is then obtained via the Perron-Stieltjes inversion formula (see 
e.g.[8,10,ll]). 

It turns out that the support of the absolutely continuous part of the measure (the 
weight function) is, for fixed A— periodic chain parameters «;„,mn, given by A, in gen- 
eral disjoint, x— intervals (bands). These are determined from the generalized Chebyshev 
polynomials of the first kind Tn{x), the so-called trace polynomials, by the condition 
|7jv(a;)| < 1. These bands coincide with the support of the spectral density (or density 
of states) of infinite chains. The weight function is, however, not given by the density of 
states. In general, a discrete point measure (Dirac 5— function measure) is also present. 
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Each gap between the N intervals which support the continuous measure may contribute 
one point (which may he on one of the band boundaries). 

In the case of Fibonacci chains the Hmit ^ oo is supposed to correspond to the 
chains based on the quasi-periodic binary Fibonacci sequence. At present this limit is 
beyond our control. 

The connection of the absolutely continuous part of the S and S measures to the 
inverse of the imaginary part of the diagonal input Green's functions for the periodic 
problem is given in an extra section. The differential spectral density (or differential 
density of states) is recovered, as usual, from the imaginary part of the average over the 
elementary A?"— unit of the diagonal (or local) Green's functions. 

In order to familiarise the reader with our notation we start with the dynamical 
equations for the displacements Qnit) = q-n exp{iu)t) at site number n 

Qn+i - 2{^-^^ - '^x)qn + Kqn-i = 0, nGZ (1.1) 

with = Kn/rrin , = Kn-i/Kn and the normalised frequency-squared x = uj'^/{2ujq). 
The spring between site number n and n + 1 has strength Kn- rUn is the mass at site 
number n. 

This recursion relation can be rewritten in terms of standard associated polynomials with 
the help of the transfer matrix method. The displacements are then given by (c/.[4a] 
(2.8), (2.10)) 

qn+i{x) = Sn{x) qi{x) - Sn-i{x) qQ{x) , neN (1.2) 
with arbitrary inputs qi{x) , q^^x) .* This expression is obtained from the transfer matrix 



■■n 



and the result 



, = , Mn-.^RnRn-i-'-Ri , (1.3a) 

qn J J 

R„ := (^<"> ) . r„(x) := 2(1^ - , (1.36) 



M„ = ( /" f-' ] , (1.4) 



where the polynomials Sn and Sn are defined by the following three-term recurrence rela- 
tion corresponding to M„ = i?„M„_i with input Mi = i?i 

= ynix)Sn-l - knSn-2 , = 0, Sq = I , (1.5a) 

= Yn+l{x)Sn-l - kn+lSn-2 , = 0, cSq = fci . (1.56) 

In the case of a mono- atomic chain (m^ = m) with equal springs (/c„ = 1) the normalised 
eigenfrequencies-squared of a finte chain with N atoms and fixed boundary conditions 



* We do not consider negative site numbers here. See ref.[4b] for the negative n case. 
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{qo = = Qn+i) are given by the zeros of Sn{x) = S'7v(2(l — x)) = Un{^ — x) , where 
U N are Chebyshev's polynomials of the second kind. These zeros are 

Xk = x[^^ = 1 - cos = 2 sin^ ^^^^ , A;= 1,2,...,7V (1.6) 

The displacements for the fc-th mode are gn+i — Sn{2{l — xj-)) Qi, for n = 1,2, ...,N, 
with arbitrary , k dependent, input qi. 

For an infinite mono-atomic chain one finds one frequency-squared band from the 
condition |T/v(l — x)\ < 1, where Tn{x) are Chebyshev's polynomials of the first kind. 
The number N of atoms in the unit cell is irrelevant because due to double zeros of 
Tn{1-x) + (-1)^+^ for X = = 1,2, A^- 1, the iV - 1 gaps degenerate. The 

x— band is, independently of A?", S = [0, 2]. In this case 5n(a^) = Sn{x) — S'n(2(l— x)), 
which are orthogonal on the inter vail [0, 2] with weight function 

w^^\x) = -^/x{2-x) . (1.7) 

Therefore there is no discrete part of the measure present in this mono-atomic case. 

The general formula for the difi^erential spectral density per particle (also called dif- 
ferential density of states) for an infinite chain with a unit of N atoms repeated pe- 
riodically is determined from the generalized Chebyshev polynomials of the first kind, 
Tn{x) := ^(SNix) — (a;)). These are the trace polynomials Mjv. (see eqn.(2.4) 

and c/.[4], eqn. (3.11)) 

for X in one of the N bands Bk, k — 1,2, N, which are determined by \Tn{x)\ < 1 and 
ordered with increasing x. Otherwise the density vanishes. 
In the mono-atomic case this becomes 

Gm{x) = G{x) = 1^4— ^ ' (1-9) 

TT ^yx{2 - X) 

for X G [0,2], and it is zero otherwise. The A?^— independence is due to the identities 
involving ordinary Chebyshev polynomials: 

1 - {Tn{1-x))^ = x{2-x){SN-i{2{l-x)f for x e [0, 2] and TV e N, as well as 

NSN-ii2{l-x)) = -Tj^'{l-x) . (1.10) 

(1.8) is, by accident, the weight function for Chebyshev's Tn{1—x) polynomials (for N = 0; 
for other N the weight is twice this). In the general case the trace polynomials {Tn{x)} 
are no longer orthogonal (see [4a], p. 5402). 

Quasiperiodic Fibonacci chains are obtained as special case with k^ = 1 , Uq/u'^ = 
-pHn) ^ -with the mass-ratio r = mi/rriQ and the quasiperiodic binary Fibonacci sequence 



h{n) = [{n + l)/(p\ - [n/(p\ , n e Nq 



(1.11) 
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In this case we use for the polynomials Sn and Sn the notation {Sn\x)} and {Sn^}. 

2. Bloch-Floquet solutions and the measure 

In this section we describe the computation of the measure with respect to which the 
polynomial system {Sn{x)} defined in (1.5a)is orthogonal, following the general method 
valid for the A^-periodic case (more generally for the asymptotically periodic case) described 
in detail in ref. [6], ch.2. The measure for the associated polynomials {Sn{x)} defined in 
(1.5b) is obtained in the same way. This method rests on the properties of the Bloch- 
Floquet solutions for the A^-periodic problem. Directly connected to these solutions is a 
map w(x) of the complex plane which enters the definition of a judiciously chosen complex 
contour integral. This integral is evaluated in two different ways i) and ii). The result 
will be the orthogonality relation and the measure can be read off. First, however, the 
necessary information on the Bloch-Floquet solutions will be given. 

The starting point is the recursion formula for the monic orthogonal polynomial sys- 
tems (indicated by a tilde) which describe N-periodic chains with Kn+N = «n and 

rUn+N = rUn- * 

S„{x) = {x- C„) Sn-l{x) - dn Sn-2{x) , 

5_i = , Sq = 1 (2.1a) 

with 

_ 1 + fen ul _ kn UjI a>^_i 

■— r, 2 ' " /I , ,2 , ,2 • \'^-^^) 

The first-associated monic polynomials {Sn{x)} satisfy (2.1a) with shifted coefficients 
Cn •= Cn+i 1 dn '■= dn+1 and the input »S_i = 0, and <So = ki. 

Orthogonality with positive definite moment functional is guaranteed for all A?^ e N by 
Favard's theorem, because dn > and the real. 

The relation between the polynomials (1.5) and the monic ones is given by (n > 1) 

Sn{x)^{-2rYl^Sn{x) , Sn{x) = {-2rkil[^Sn{x) . (2.2) 

For the Fibonacci chains with Kn = k one uses the N— periodic binary sequence h^^'' (n) 
obtained by repetition of the first N entries of the quasiperiodic sequence {h{n)} of (1.11). 
The original polynomials {Sn\x)} and {Sn\x)} for quasiperiodic Fibonacci chains, are 
then obtained in the limit N oo, keeping always n < N. 

In the A^— periodic case the transfer matrix Mn of the elementary unit satisfies: Det Mn = 
kq/ Kn — +1- From this and Mn+N = MnM^ one derives, using (1.4) 

Sn+2N{x) = 2 Tn{x) Sn+Nix) - Sn{x) , (2.3a) 
Sn+2N-l{x) = 2 Tn{x) Sn+N-l{x) - Sn-l{x) , (2.36) 



* All quantities depend on the chosen period N, uP' = (a;o, o^i, a;^_i) and k 
{k\i k2, kN) with kq = kn- In the sequel this dependence will be suppressed. 
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with the trace polynomials 

Tn{x) := ^tr Mn = ^{Sn{x) - SnM^)) ■ (2-4) 

In the Fibonacci case we use the A^— periodic binary sequence {/iAr(?i)} obtained from 
(1.11) by taking hN{n) — h{n) for n = 1, 2, and defining hN{n + N) — hN{n). In this 
case the polynomials Sn, Sn-2, hence T/v, coincide with the corresponding quasiperiodic 
ones. 

(2.3) implies for the general solution for the displacements (1.2) 

qn+2N+i{x) = 2 Tn{x) qn+N+i{x) - qn+i{x) . (2.5) 
The two independent Bloch-Floquet solutions are, for arbitrary input qo{x) and qi{x), 

qn,±ix) := qn+N{x) - Xn,±{x) qn{x) , (2.6a) 

with 

XnA^) ■■= ^n{x) ± ^J{TN{x)f-l , (2.66) 

which are the two eigenvalues of Mjv- With this definition the periodicity condition (2.5) 
yields qn+N,±{.x) = Xn,^:{x) qn,±ix) ■ Xn^z^{x) = exp ij3N,^{x) for x values in any of the 
N intervals (bands) Sfc , k= 1,2, ...,N, defined by * 

\rN{x)\<l , (2.7) 

which is necessary for harmonic vibrations. Therefore, 

qn,±{x) ^ {XN,^{x))'^ifn{x) , n G Nq , (2.8) 

with periodic <^„_|_7v(a^) = '^n{x) , proving that (2.6a) defines the Bloch-Floquet solutions. 
The integrated spectral density (or integrated density of states) is a continuous function 
given by appropriate branches of the Bloch-Floquet phase (3n{x) — cos~^Tn{x) for band 
values X and constant pieces in the N — 1 gaps between the bands. The differential spectral 
density (or differential density of states) (1.7) follows from this after differentiation. 
In order to compute the measure with respect to which the polynomials {S„,{x)} are 
orthogonal one considers (see ref. [6]) the following map for a; e C 

w{x) ^ {Xn{x))^ , (2.9) 

with the sign of the square-root in (2.6b) chosen such that for real x G B^, k = 1, N, 
Xn{x) runs along the unit circle from -|-1 to —1 for odd numbered bands and from —1 to 

* From (1.5) one finds 3^(0) — I+kq Yl7=i l/^j ' '^^(0) — ki+Ko Yl^=2 V'^i ■ Therefore, 
in the N— periodic case Tn{0) = 1 holds for all chains (kq = kn)- For the proof of the 
reality of the zeros of Tat (a;) ^ 1 see ref. [6], lemma 2.2. 
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+1 for even numbered ones. This requires the sign (—1)'^+^ for x e Bk- Xn{x) is real for 
X in the gaps Gk, /c = l,...,A^— 1, and we choose the sign of the square- root as ( — 1)'^ such 
that |AAr(x)| > 1 holds, x G {xmax, oo), with the maximal band value Xmax^ is taken as 
gap Gn with the sign choice (—1)^. Such a sign choice poduces for u!{x) of (2.9) a path 
like the one shown for the Fibonacci chain example A/" = 5, r = 2 in fig. 2, if x runs along 
the real axis from to Xmax- Finally, in the z = l/w{x) plane a closed contour is obtained 
if X runs also backwards from Xjnax to x = 0. This contour is shown for the example in 
fig. 3. The starting point is at 2 = +1. In figs. 2 (resp. 3) the S^, (resp. Bf^) and G;^ (resp. 
Gf^) labels indicate the images of the bands B^ and gaps Gk under the map w(a;) of (2.9) 
(resp. z = l/w(a;)). That the map w is of importance is clear from the Bloch-Floquet 
solution (2.8). 

For the computation of the measure for the polynomials {Sn{x)} it was found in ref.[6] 
that one should consider the plane contour integral for m, n = 0, 1, 2, ... 

with g the inverse map to z{x) = l/w{x), viz x = g{z), and the Bloch-Floquet solution 
Qn+i,±ix) defined by (2.6a) and (1.2) with the choice qi = 1 and go = 0. The sign is chosen 
as described above in the definition of w{x). In addition, for Xmax < x < +00 the sign 
choice is {—)^- The contour V (with negative orientation) is the union V — V ^" U Vq" 
with 

F^// := {z\z = expid , 9 ±^k , k = 1,2, N - l} (2.12) 

and Tq" are the 2{N — 1) closed curves around the images of the gaps Gk- For 
the Fibonacci chain example with = 5, r = 2 see figs. 1 and 3. For later use note that 
because gn+i,±(^) satisfies the same recursion relation as Sn{x) one finds, after comparison 
of the n = — 1 and n = inputs on both sides, (remember go = 0) 

qn+i,±{xi) / qi,±{Xi) = Sn{Xi) , (2.13) 

with the zeros Xj, i = 1,2, N — 1, of S^-i. The contour integral is computed in two 

ways i) adding the F^" and Tq" contribution, both transformed to the x— plane, and ii) 

evaluating the 2;— plane contour integral with the help of the residue theorem. 

In the first step of i) one computes the Tb" contribution using, for x & Bj, with sign choice 

(-1)^+1, 

qn+l,±{x) = Qn+lM^) + {'^/>^N,±{x) - Xn,±{x)) Sn{x) , (2.14) 

following from (2.6a), and with the choice gi = 1, go = in (1.2) one has g„+i = «S„. After 
some rewriting one finds (see appendix Al for details) 

Ib = -—I iy^^i^)-l^(^))^rn{x)S^ix) 

2m Jr „ Sn-i{x) 

B 
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with X = g{z), and B'^ denotes the lower half of the punctured unit circle in the 2;— plane. 
This can be restated in the x- plane (see fig.l for the Fibonacci case N = 5,r = 2): 

1 f Smjx) >S.(x)(-l)^+Vl - iT^{x)Y _ ^2.16) 



Jb Sn-i{x) 

The second step of the computation is that along Tq"- This integral is again evaluated 
in x-space, such that the — 1 (positively oriented) closed contours around the gaps Gk 
are picked up (c/. fig.l for the = 5, r = 2 Fibonacci case). Sn-i{x) has a simple zero in 
each of the — 1 gaps. This follows from orthogonality and a lemma (Ya.L. Geronimus, 
M. Kac and P. van Moerbeke, Lemma 2.2, page 46 in ref.[6]) about the interlacing of the 
zeros of Tjv T I5 which are real, and <Sjv-i (and <SAr_i). Note, that <Sjv-i zeros may occur 
at one of the band boundaries, coinciding then with zeros of TJv — 1 or T/v = 1. The residue 
theorem now yields with the help of (2.13) 

/N-l , N 

dx{ J2 ^{^ - Sm{x) Sn{x) . (2.17) 

k=l '^AT-ll^j 

with the zeros Xk = ^^'^ of S^-i. 

Computation ii) of (2.11) is done in the 2— plane using the residue theorem. The only 
possible pole inside F, coming from (/'(-z), occurs for a; — > 00, i.e z = 0. The zeros of Sn-i 
are outside of F. Therefore only the large x behaviour of the integrand is of interest. 
For large x the <S— polynomials behave like (see (2.2)) 

S^{x)r.{-2rl['^x^ . (2.18) 

i=i ^« 

The asymptotics of qn+i,(-)if{x) cannot be found from its definition (2.6a) directly. Fol- 
lowing ref. [6], it is found from the fomula, based on (2.6a) and (2.5) 

qn+l,+ {x) qn+l-{x) = {qn+l+N{x)f - qn+2N+l{x) Qn+lix) ■ (2.19) 

The r.h.s. can be rewritten, using the m— th associated polynomials sji"^ defined by 

5^)(a;) = r„+^(x)5(!!U^)-^n+m5S(^) , (2-20) 

with inputs S^"^ = and iSq™''* = YliLi m G N and Sq^^ — Sq = 1. Because 

satisfies the recursion relation oiq^^i, with the two independent solutions gn+i and 
gn+i+AT, with Wronskian g„+i g„+jv - Qn+i+N Qn satisfying 

n 

W{qn+i,qn+i+N)/Ylki = W{qi,qN+i) = qiqN - QaQN+i , (2.21) 



i=l 
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one finds (for general qi and qo) 

'^i-Ll^) = yx./ ^{im+Nix) qn+i{x) - qm{x) qn+i+N{x)} . (2.22) 

^ [Qi, Qn+i) 

Letting n ^ n + N and m — > n + 1 , one obtains the r.h.s. of (2.19): 

qn+i,+{x) qn+i,-{x) ^W{qi,qN+i) S'^:!:p{x) . (2.23) 

The asymptotic form of qn^i^(^_^N can now be inferred from the one of g^_|_i^^(_)]v which 
follows without difficulty from (2.6a) *. For gi = 1, go = one finds that the leading term 
for large x is 

/ \ n+l n+1 2 

«=±HzHf) ^ Yl fc,/((-2)»+' n ^ ■ (2-24) 

The residue theorem can now be applied to the contour integral (2.11) in the 2— plane. 
g'{z) ~ —C I for small z due to (7(2) = x and ^ = Ijwix) ~ C/x for large x. The 
capacity C drops out in the calculation of the residue for 2; = 0. The result is 

-^m,n = +2 ^Jl W ^n,m • (2.25) 



n- 



Combining both ways of calculation i) and ii) one ends up with the normalized measure 
da for the orthonorma/ polynomials with sq = 1 



sn{x) := (-1)" J f,' }^ Sr,{x) , (2.26) 



where the factor (—1)" has been inserted to guarantee positive leading coefficient. 

J Sm{x) Sn{x) da{x) = 5m,n , m, n G Nq (2.27) 



where 



(2.28a) 



fe=i 



^ 1 2 ^ (-iryr^M , . e , . = 1,2, iV, (2.286) 

TTKiUjf Sn-i{x) 

qi,(-)>^ix) = Snix) - XM,(-)>^ix) , k = l,2,...,N-l, (2.28c) 

XN,i-)4x) =Tn{x) + ^/iTNix))^ , xeGk, (2.28d) 



* The leading coefiicient of S^^^^ (x) is, for n e Nq, (-2)^-^ nSi nr=n+2 ^o/^f ■ 
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The absolutely continuous part of the measure, w{x) dx, vanishes outside the N bands 
B = UBk- It is non-negative because the sign of 5iv-i in band is ( — 1)'^"'"^, due to the 
fact that iS7v-i(0) > +2 and the interlacing property of its zeros with the one of Tjv =F 1- 
The Dirac-measure lives on the A'" — 1 zeros of Sn-i- The fact that also this measure 

is non-negative will be discussed in the section 4. 

A similar computation can be performed in order to find the measure for the associated 
orthogonal polynomial system {<S„} of (2.2b). One puts in (2.1) qi = and qq = —1. In the 
integral (2.11) one uses S instead of S and replaces qn+i,± by qn+2,± '■= Sn+N — ^n,± S^. 
In place of (2.13) one uses here qn+2,±{,Xi) / q2,±{xi) = Sn{xi)/ki with the zeros Xi of 
Sn-i- The normalized measure for the orthonorma/ polynomials with positive leading 
coefficient and So = 1 



kl V ^n+2 n, = l kj 



is then 



d»(.) = dx - ^'^(x-lf-')) 24|l ^^.y'Wf ' , (2.30a) 

7ra;|/c2 Sn-i{x) 

q2,{-)k{x) ^ Sn{x) - ki XN^^-)k{x) , xeGk, /c = l,2,...,iV- 1, (2.30c) 

where A7v,(_)te is given in (2.28d), and are the zeros of Sn-i{x). 



3. Continued fraction, Stieltjes inversion formula, and the measure 

The second method to compute the measure for orthogonal polynomials with periodic 
recursion formula coefficients relies on the fact that the continued fraction accompanying 
this recursion formula is the Stieltjes transform [8] of the orthogonality measure [lb, 9]. In 
the periodic case the continued fraction can be given explicitly, and the measure is then 
determined with the help of the Perron-Stieltjes inversion formula [10,11]. 
The continued fraction which belongs to the recursion formulae for the (not necessarily A^- 
periodic) associated orthogonal polynomials {Sn{x)} and {Sn{x)} (see (1.5a) and (1.5b)) 
is 

kiLuf ^ . ki \ k2 \ kn \ /o i\ 

X[x) = I / X I / X - I / X , (3.1) 



2a;2 ' I Yi{x) \ Y2{x) \ Yn{x) 
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where Yn{x)) is given by (1.3b). The factor —kiUi/2LiiQ has been introduced for later 
convenience. The n-th approximation to this continued fraction is for n > 1 * 

h^^f /X ki \ k2 \ kn \ <Sn-i{x) , . 



2u^o ^ I Yi(a:) | Y2{x) - | r^(x) 5,(x) ' 

which follows by induction, using the recursion formulae. If the continued fraction con- 
verges x(x) := lim^^oo Xn{x) . A fundamental theorem (see e.g. [9], theorem 2.4, or [lb], 
p. 90) states that x{x) is the Stieltjes transform of the measure 

f+°° da(t) ^ , , ^ , , 

X{x) — / X suppyda) . (3.3) 

J —oo ^ 

Here da is the real, positive, and normalized measure for the orthogonal S— polynomials 
(1.5a) (c/. [9], eqs.(2.1) to (2.5)) 

Sn{x) Sm{x) da{x) ^ -^^W kj 5n,m = —Sn,m , m, n G Nq (3.4) 

The Perron-Stieltjes inversion formula (see e.g. [9, 10]) for a real measure of bounded vari- 
ation is 

1 f*^ 

cr{t2) - cr{ti) = lim / Im x{t + iv) dt , (3.5) 

with cr(tfc) := ^{a{tk + 0) + cr{tk — 0)) for A; = 1,2. xix) = xi^) for x e C, and x is analytic 
for non-real x. 

In the N— periodic case xix) can be calculated explicitly as follows (c/. [7]). Consider, for 
fixed X, N and parameters /c„, o;^, the map in the complex ^— plane 



ki I /c: 



JNix; z)=z' ^ ' - ^ , (3.6) 

with l^(a;) given in (1.3b). By induction, with the help of the recursion formulae, one 
finds 

Jn{x;z) = — -r- ^ . . , (3.7) 

<Siv-i(a:) z - Sn{x) 



The n-th approximation to the continued J-fraction is for n > 1 (see (2.1) and (2.2)) 

ki \ d2 \ dn \ , . <Sn-i{x) 

1 i ••• ~ 1 = Xn{x) = — • 

\X-Ci \X-C2 \X-Cn Sn{x) 
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which is for real x a SL{2,Ti) Mobius transformation due to 

1 = Det Mn{x) = -Sn{x) Sn-2{x) + Sn-i{x) Sn-i{x) . (3.8) 

Because of A^— periodicity —kiU!fx{x)/'2i^o is a fixed point of the map (3.6), and can be 
computed from (3.7) as 

-^XN,±{^) = {l{SNix)+SM-2{x)) T VCTNixW - l}/5iv-i(x) , (3.9) 

where (3.8) was used, and Tn{x) is given in (2.4). For |77v(a;)| < 1, which determines N 
bands Sfe, A; = 1, 2, N, Xn,±{x) becomes complex. For the gaps between the N bands, 
Gki A; = 1, 2, A^— 1, it is real. See fig.l for the Fibonacci case N = 5,r = 2. Introducing 
^N,±ix) given in (2.6b), this is rewritten as 

k up' 

~":r^Xiv,±(a^) = {Sn{x) - Ajv,±(a;)}/5jv_i(a;) . (3.10) 

Xiv,(-)^(^) is for large x proportional to 1/a;.* The sign choice for x values in the bands 
and gaps will be given below. 

The measure can now be determined from the inversion formula (3.5). The absolutely 
continuous part of the spectral measure, daacix) = w{x) dx, is found from 

d 1 

w(x) = -;-(Jacix) — lim Im xix + ir]) < oo . (3-11) 

dx TT »7->+0 

x{x) = Xn{x) has to be defined such that x{x) — x{^) holds for complex x. This single- 
valued function is called xix)- w{x) lives therefore on the bands Bk and coincides with 
(2.28b), after the sign of the square-root {i.e. of the Riemann sheet) has been chosen such 
that w{x) is non- negative. See section 4 for the proof that the sign choice in (2.28b) leads 
to positive w. 

The discrete part of the spectral measure (sum of Dirac 5-functions) originates from the 
simple poles of x{^) with their residues determining the height of the jumps of the measure. 



dao....(x) = ir'') ^"'r! " tf^ ■ (3.12) 



where the sign of the square-root in \n{x) of (2.6b) has been chosen as (— l)'^ for x e G^, 
like in the calculation leading to (2.28c). This choice will be seen in the next section to 
produce positive jumps in (7{x) at the zeros of <Sjv-i . For finite N and parameters 

{cu^} there are no other singularities of x{x) which could contribute to the singular 
part of the measure (c/. [12], ch. XIII, p. 140). 



* The a;— asymptotics cannot be found from (3.10). One uses 
(<Sjv — Ajv,(-)Jv)/<SAr-i = Sn-i/{Sn — Ajv,(-)^+i) which is identity (2.23) with = , 
qi = 1, (2.6a) and (1.2). 
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The measure for the orthogonal {Sn{x)} polynomials can be calculated in a similar fashion 

(2,) 

by taking into account also their first associated polynomials {Sn (x)}, defined in (2.20). 
Note that because of the cyclic property of the trace and A?"— periodicity on has: 

fN{x) := {Sn{x) - sf_^{x))/2 = ^tr Mn = hTNix) (3.13) 

whith Mn '■= Rn+iRn ■ ■ ■ R2 ■ More details are found in appendix A. 2. 
We have thus reproduced the results of the previous section. 

4. General remarks and Fibonacci chain examples 

We first state a simple conclusion concerning the discrete measure (3.12), which will 
show its non-negativity. With the mentioned sign choice the numerator of the k-th term 
can be rewritten, with (2.6b), (2.4), and (3.8), like 



Sn — — -^{Sn + SN-2) + {~^)'^~^^'^ {2^^^ — Sn-iSn-i • 



(4.1) 



Evaluated at the zero Xk = of »Sjv-i this becomes 



(«Siv-Aiv)| ^ i^{-i^N + SN-2)\xk) {Sn + SN-2)\xk for k even 
&{{Sn + SN-2)\xk) i^N + SN-2)\xk fork odd 



(4.2) 



with the step function 0(x). Due to (3.8) {Sn + SN-2)\xk — {{^N{xk))'^ — ^) /SN{xk)- The 
final result for the discrete measure is 

■i-KcW = v^T. <s(-«r-") e((-)'+'(5«(x)+5„_. (.))) ■s.vW + ^.v-.M ^ _ 

(4.3) 

which is always non-negative because the signum of «5jy^_i(^^'^ is (— I)*'- 
Two remarks are in order. 

i) There will be no contribution to the discrete measure from those zeros ofSN-i which sat- 
isfy Tn{^1!^ ^'*) = ±1 .In this event the zero oiSN-i coincides with one of the band bound- 
aries, and from the definition of Tn and (3.8) one finds SN{^i^~^^)— —SN-2{Ck^~^^) ~ 

ii) A comment on band degeneracy. A gap Gk will disappear whenever Tn + (— 1)^"'"''^ 
has a double zero at, say, Xk. Because the zero of Sn-i lies in the gap or on 
one of the adjacent band boundaries one has Xk — ^^'^"^^ For the same reason the k- 
th zero of »Sjv-i is then also Xk- Due to (3.8) SN-2{xk) = — l/<5Ar(xfc), and therefore 
SN{xk) + l/SN{xk) = 2{-l)'' (definition of Tn). Thus SN{xk) = (-1)'' = -SN-2{xk) , 
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and there will be no contribution to the discrete part (4.3) of the measure from such 
disappearing gaps Gk- * 

Next, we consider some examples of Fibonacci chains. 

In the AT— periodic Fibonacci case fcn = 1 and Yn{x) = 2(1— r'^^^^^a;), with the AT— periodic 
binary sequence {hN{n)} obtained by continuing the first N entries of {h{n)} given by 
(1.11) periodically. All polynomials will now depend on N and the mass-ratio r = m^i/rriQ. 
We shall use non-script symbols for these polynomials and the r labels will be clear 
from the context. 

a) First we check the mono-atomic case r — 1. The results have already been given in 
the introduction. Remark ii) applies for all A?" — 1 disappearing gaps. 

b) Put N = 2, r = niA/mB 7^ 1, i-e. AB-chains. We first show that there is no discrete 
part (4.3) of the measure. Here remark i) applies. The zero of Si{x)= 2(1 — rx) (see [4a]) 
is xi = 1/r, and because 

T2{x) = 2rx^ - 2(1 + r)x + 1 , S2ix) = 4:rx^ - 4:{1 + r)x + 3 , (4.4) 

one has S'2(l/r) = T2{l/r)=—l=—So. ** The weight function is given by 



2r / 1 11 

w{x) = — \ —x{x )(a;— l)(a;— 1 ) /\x 1 , (4-5) 

TT \ r r r 

for X in any of the two bands 

r>l : Si = [0,-] , S2 = [l,l + -] , (4.6a) 
r r 

r<l : Si = [0,1] , S2 = [-,l+-] . (4.66) 

r r 

The weight function (2.30) for the associated polynomials {Sn{x)} is found to be 

w = -1^—4 w(x) , (4.7) 
r \x — 1\ 

for x in the bands (4.6) and zero otherwise. 

c) Put = 3, r = mA/mB 7^ 1, i-e. infinite ABA-chains. Now the discrete measure is 
present because zeros of S2 satisfy 2rx^: = 1 + r ^ ^/(r — 1)^ -|- r and S3{xzf)+Si{x^) = 
2{r - l)x^ ^ 0. Therefore 

J _ 4r J (r — 1) a;_ (5(a; — a;_) da; for r > 1 . 

~ ^ (r- l)2 + r I (1 - ^) x+ -x+)dx for < r < 1 ' ^ 



* In the A^— periodic Fibonacci case an example is N=:6 ( [ABA)'^ chains) where one 
finds double zeros of + 1 at the zeros of for /c = 1, 3, 5. 

** This is an example where a zero of Sn-i coincides with a band boundary without 
having band degeneracy (no double zero of Tjv T !)• It seems to be a counterexample to 
one part of the statement found in [6], lemma 2.2, top of page 47 (the "if part). 
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The bands are, depending on sign{l — r), 

r > 1 : = [0, -] , B, = [-6_(r), ^] , ^3 = , (4.9a) 

r<l: B, = [0, , i?2 = ^] , S3 = [|:, ^6+(r)] , (4.96) 



with 6±(r) = r + I ± y — r + |. 

The weight function (2.28b) becomes ior x e Bi^ , k = 1, 2, 3, 



2 , V-^(^-^)(^-|:)(^-i)(^- j7^+(r))(x-^6+(r)) 
^^^^ - TT^ (-l)^+i(x - id-(r))(x - ici^r)) ' ^'-'^^ 

with (i±(r) = 1 + r ± v^r^ — r + 1. 

Put = 4, r = m^/m^ 7^ 1, z.e. ABAA-chains. There is no contribution to the discrete 
measure because one finds (see [4a]) for the zeros of S^, vizxi = 1/r, 2rx± = l+r±-\/p"+T, 
S4{xk)+S2{xk)=0 for k = 1,±. Like in the N = 3 case one can give the bands and the 
weight function exphcitly. 



5. Green's functions and the measure* 



The purpose of this section is twofold: i) to define the Green's functions of the 
general A?"— periodic problem which satisfy specific boundary conditions appropriate to the 
use of the Bloch-Floquet solutions encountered in section 2. ii) to find the orthogonality 
measures computed in this work from these Green's functions. We shall also give the 
relation of the differential spectral density (or differential density of states) (1.8) to theses 
Green's functions which turns out to be the standard one. This should clarify the difference 
between the absolutely continuous part of the measure and the spectral density. 
The Green's functions Gn,m{x) for the A?"— periodic problem are defined for n,m & Nq by 



Yn+l{x) Qn,m{x) - Qn+l,m{x) - kn+1 Qn-l,m{x) = Sn,m ■ (5.1) 

The F— coefficient is defined in (1.3b). The inputs are Q-i,m{x) and Qo^mix). The Green's 
functions with the proper boundary conditions turn out to be the ones constructed from 
the Bloch-Floquet solutions (2.6). This solution is of the type 

Gn,m{x) = 0"m{x) Qmax{Ti,m) + l, {-!)'' {^) 9min(n, m)+l, (-1)^+1 (^) . (5.2) 

* This section is inspired by a paragraph found in ref.[7] for periodic problems. 
** The dependence on A^, {kn}, {^^n} suppressed. For the monic polynomials (2.1) 
one uses {x - c^+i) Gn,m{x) - dn+i Gn-i,m{x) - Gn+i,m{x) = Sn,m ■ The relation 

(nr=lSt)^n,m = i-'^r-'^+'intV Sn,m holds. 

i i 
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The sign choice pertains to gap Gk, for which qn+\^{^-i)k — > for a; e Gk and n ^ oo 
(sec (2.8)). The coefficient is found from (5.1) putting n = m and using the recursion 
formula for the qn,±- For general input qi and go one finds with the Wronskian (2.21) 



m+l ^ 



ki Gn,m{x) — 



W{qi,qN+i) (A7v,(-i)fc(2;) - A^^(_i)fe+i(x)) 

■ Qmax{n,m) + l,{-l)'^{x) ?mm(n,m) + l,(-l)*+i (^) • (^-3) 

This is the Bloch-Floquet Green's function vanishing for a;— values in gaps for n — > oo with 

fixed m and vice versa. 

Our interest is in the diagonal Green's functions Qn,n{x) which are in fact go and gi 
independent. This is due to identity (2.23) which shows that the Wronskian drops out. 
The final result can be written in terms of the m— th associated polynomials {S^\x)} 
defined in (2.20) like 

n+l . 

(H GnA^) = S^^^l\x)/(2^J{T^{x)f - l) . (5.4) 

i=l 

The sign of the square root depends on the gap and band number. For the A;— th gap, G^, 
it is (— l)'^, for the A;— th band, it is (—1)'^+^ in accordance with the remarks found in 
section 2. In particular, one has for the input quantities 

g-i,-i{x) = ^Sn-i{x)/^{Tn{x)Y -1 (5.5a) 
GoA^) = ^^Sn-i{x)/^ {rN{x)f -1 . (5.56) 

The imaginary part of these input Green's functions are inversely related to the weight 
functions computed in this paper. For x&Bk, k = l,2,...,N, one finds 

Img-i,-i{x) := lim /m 6?-i,-i(a; + ir/) = 1{-1)''Sn-i{x)/ Jl - (rjv(a;))' . (5.6) 

In the gaps the imaginary part is zero. This shows that the absolutely continuous part 
of the {Sn{x)} measure (the weight function), which lives on the bands, is essentially the 
negative inverse of the imaginary part of the _i Green's function. 

-/m^_i_i(a;) = u}Q/{7ru)lw{x)) , (5.7) 

with (2.28b) and x in the bands. Similarly 

-IrnQoAx) = u>l/{nu>lk2w{x)) , (5.8) 

with (2.30b). The average over the elementary A?"— unit of the chain for the diagonal 
Green's functions gn,nix), belonging to the monic polynomials {Sn{x)}, can be computed 
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with the help of the Christoffel-Darboiix identities for the Bloch-Floquet solutions (2.6). 
These identies follow from the recursion formula (c/'.refs.[l]), and they are (remember that 
kik2 ■ ■ ■ kjsf = 1 in the A^— periodic case) 

N-l 2 

n=0 lli=l '^i 

Qn,+ {x) q'N+i-{x) - qN+i,+ {x) q'N-{x) + qi,+{x) qo_{x) - qo,+{x) q[_{x). (5.9) 

There is an alternative version of this identity where the derivative acts on the left factor 
and an overall minus sign appears. This happens because the confluent Christoffel-Darboux 
identity implies 

qN,+ {x) qN-{x) - qN,+ {x) qN+i-{x) - qi,+ {x) qo,-{x) + qo,+{x) qi-{x) = , (5.10) 

which is in fact the Wronskian identity analogous to (2.21) for W{qN+i,+ ,qN+i,-)- The 
connection between Qn,nix), for the monic case, and Qn,nix) is 

GnA^) = -2 ^ gn,n{x) . (5.11) 

Using (5.4), (2.23) for the (n + 1)— th associated polynomials, identity (5.9), and the 
definitions (2.6) with (1.2), a lenghty calculation shows that for general inputs qo{x) and 

qi{x) 

N-l 

E ^nA^) = r^{x)/^l{r^{x)f - 1 , (5.12) 

n=0 

with the sign convention for gaps and bands mentioned earlier. Therefore, the imaginary 
part produces the differential spectral density (or differential density of states) known from 
the differentiation of the Bloch-Floquet phase: 

1 1 _ 

- lim Im — Y] Gn,n{^ + '^V) = Qn{x) , (5.13) 

TT ri^O+ IS ^ 

n=0 

given by (1.8) for x values in the bands and zero otherwise. This result corroborates the 
choice of the Bloch-Floquet Green's functions. As a by-product we find from (5. 12), (5. 11), 
and (5.4) the identity 

iV-l 2 

E 2 nVi, = -^iv'(-) , (5.14) 

for the A?"— periodic case of the general associated polynomials (2.20) which collapses to 
the well-known identity (1.10) for Chebyshev polynomials for the mono-atomic case. 
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Appendix 

Al. Derivation of eqn.(2. 14) (c/. Ref.[6]) 

In the 2;— plane Ib is given by (2.11) with F = Tb" defined by (2.12). The sign 
choice of qn+i,± depends on the band Bk, k = 1,2,...,N, where it is (—1)'^+^. qq = 
and (/i = 1 in (2.6a) with (1.2). The relation (2.14) with the appropriate sign choice for 
Xn{x) = A7v,(-i)fe+i (2:) is used in order to rewrite g^+i.i of (2.11). The piece with g„_|_i^rp 
is, after a change of variable, seen to be the negative of the original integral Ib with the 
relevant qn+i,± choice for x e -Bfc- The same change of variable is used to rewrite the 
second piece coming from (2.14) as twice the integral over only half of the contour, namely 
over T^" in the lower half of the 2;— plane. 

A2. Continued fractions for the {<Sn} measure calculation 

The {Sn{x)} recursion relation is given by (1.5b). Their first associated polynomials 

(2) 

are {Sn (x)} defined by (2.20). The normalised measure da obeys 

n+2 ^2 ^ 
Snix) Sm{x) da{x) ^ -^Y\ ~'^^ri,m = — dn,m , m, n G Nq . {A.l) 

The continued fraction X(x) which is related to this measure like in (3.3) has approximants 

X (x) ^"-^^^^ ^"^^ ( ' ' ^-+1 I ^ M2) 

^ " Sl{x) " k2ul\\ Y^{x) I ^3(0.) - \Y^^,{x) ■■■) ■ ^""-'^ 

The tilde quantities are monic polynomials. Sn^ = Sn The corresponding Mobius trans- 
formation is 

J z)^z' = ^ - ^3 I - - ^iv+i I ^ ^N^-2{^) ^ - ^N^-l{^) 

^ 1^2(3;) 1^3(0:) \Yn+i{x)-z Sn-i{x) z - Sn{x) 

(^.3) 

The fixed point solution can be written like 
7 2 

^'^''-Xi,^^{x) = {SM{x)-k^XN,±{^)]/SN-i{x) , (A4) 



2a;2 



where Xn,± '■= T'Nix) ± y {T^Y — kf — kiXN^± was used which follows from (3.13) and 

(2.6b). The measure da is then computed like in (3.5) from X. With the definition (2.29) 
of the {sn} polynomials one finds / da Sn{x)sm{x) = 5n,m with da given in (2.30). 
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